HOLOMORPHIC SPINOR OBSERVABLES 
IN THE CRITICAL ISING MODEL 



DMITRY CHELKAK^'*^ AND KONSTANTIN IZYUROV^'^ 



Abstract. We introduce a new version of discrete holomorphic obscrvables 
for the critical planar Ising model. These observables are holomorphic spinors 
defined on double covers of the original multiply connected domain. We com- 
pute their scaling limits, and show their relation to the ratios of spin correla- 
tions, thus providing a rigorous proof to a number of formulae for those ratios 
predicted by CFT arguments. 



1. Introduction 

1.1. The critical Ising model and Smirnov's holomorphic observables. 

The two-dimensional Ising model is one of the most well-studied models in statistical 
mechanics. Given a discrete planar domain V,^ (a bounded subset of the square 
grid) , the Ising model in can be viewed either as a random assignment of spins 
to the faces of , or a random collection of edges of fi"^, with an edge drawn 
between each pair of faces having different spins. The partition function of the 
model is given by 

^ ^ TifMn] or 
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respectively, where J-{^^) denotes the set of faces of il^ and Conf(il'') is the set 
of subgraphs S of such that all vertices of have even degrees in S. We 
refer the reader to Section [2] for a more detailed discussion and notation. We 
will be interested in the properties of the model at the critical temperature T = 
21og~^(V2 + 1), which corresponds to x = \/2 — 1. This value of x will be fixed 
throughout the paper. 

Discrete holomorphic observables, also called holomorphic fermions or fermionic 
observables, were proposed by Smirnov in |Smi06j as a tool to study the critical Ising 
model, although similar objects appeared earlier in |KC71| and |Mer01j without 
discussing corresponding boundary value problems. Since then, these observables 
proved to be very useful for a rigorous analysis of the planar Ising model at criticality 
in the scaling limit when fl^ approximates some continuous domain ft as the lattice 
mesh S tends to zero. 
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Recall that Smirnov's fermionic observable is defined as 



(1.1) 



where Confa^z is the set of edge subsets S, such that S can be decomposed into a 
disjoint collection of loops and a simple lattice path 7 connecting a boundary edge 
a to the midpoint z = Ze of an inner edge e; w(7) is the winding number of 7; and 
na & {±l,±i} denotes the orientation of the outgoing boundary edge a. With this 
definition, the observable has been shown to be discrete holomorphic and satisfy 
Riemann-type boundary conditions 



This led to a proof of its convergence to a conformally covariant scaling limit |CS12j . 

This result has been the main ingredient of the recent progress in rigorous under- 
standing of conformal invariance in the critical two-dimensional Ising model. The 
martingale property of F{a,z) (see further details in |US12j ) allows one to prove 
convergence of the Ising interfaces to the chordal Schramm's SLE3 curves. Using a 
slightly different version of this observable, Hongler and Smirnov |HS10j were able 
to compute the scaling limit of the energy density in the critical Ising model on 
the square grid, including the lattice dependent constant before the conformally 
covariant factor. This result was later extended to all correlations of the energy 
density field and certain boundary spin correlations [HonlOj . 

At the same time, similar observables proved to be very useful in the analy- 
sis of the random cluster (Fortuin-Kasteleyn) representation of the critical Ising 
model |Smi06l iRCOHl ISmilOi [CST2I IDHNl Ij . In particular, it was shown by Beffara 
and Duminil-Copin jBDlOj that they can be used to give a short proof of criticality 
of the Ising model at the self-dual point. 

Many of these results generalize beyond the case of square grid approximations. 
Thus, convergence of fermionic observables has been proven for isoradial lattices 
[CS12 ■ which reappeared in the connection with the critical Ising model in the 
paper of Mercat [MerOlj . This proved the universality phenomenon, i.e., the fact 
that a microscopic structure of the lattice does not affect macroscopic properties 
of the scaling limit. Moreover, discrete complex analysis technique developed in 
[CSllj and |CS12j provides a general framework for such universal proofs. 

On the other hand, one of the most natural questions about the Ising model - 
the rigorous proof of conformal covariance of spin correlations in the scaling limit 
— remained out of reach until recently. The goal of the present work is to introduce 
a new tool - spinor holomorphic observables ~ that allows to attack this problem. 
In particular, we prove convergence of ratios of spin correlations corresponding to 
different boundary conditions to conformally invariant limits. In a subsequent joint 
paper with Clement Hongler |CHI12| , using a more elaborate version of the spinor 
observables, we prove conformal covariance of spin correlations themselves. 

1.2. Spinor holomorphic observables and ratios of spin correlations. In 

this paper we extend the study of fermionic observables to the case of multiply 
connected domains. Given a double cover zu : fl^ il^ of such a domain, we define 
the observable F^{a, •) : 17"^ — C by 





(1.2) 



SeConf, 
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where a, z £ Q , but the sum is taken over the same set of configurations as before; 
l{s) is the number of loops in S that do not hft as closed loops to il^ , and 
is the indicator of the event that 7 lifts to as a path running from a to z (and 
not to the other sheet), see Section [s] for detailed discussion. In other words, we 
plug into ( |1.1[ ) an additional sign that depends on homology class of S modulo two. 
It is worth to mention that our observables should be closely related to the vector 
bundle Laplacian technique applied to uniform spanning trees and double dimers 
by Kenyon |KenlO[ IKenll] , although at the moment we do not know any exact 
correspondence of that sort. 

Our main observation is that F^{a,z) are discrete holomorphic and satisfy the 



boundary conditions (1.2 1, just like Smirnov's observable F(a,z). The definition 
implies that F^(a^ z) = —F^{a,z*), \i z ^ z* belong to a fiber of the same point; 
hence, we call F^ holomorphic spinors. 

To describe the scaling limits of F^(a, •), we will introduce the continuous holo- 
morphic spinors fm{a, •). Roughly speaking, th ese a re fundamental solutions to the 



continuous Riemann boundary value problem (1.2) on the double-cover O, with a 
singularity at a and the property /ct(-z) = ^ft^{z*). Postponing precise definitions 
until Section [sj we now state our first main result (see Theorem 3.13): 



Theorem A. Suppose that ft^ is a sequence of discrete domains of mesh size 
5 approximating (in the sense of Caratheodory) a continuous finitely connected 
domain Q, and that a* € dil^ tends to a G dil, as S ^ 0. Then there is a sequence 
of normalizing factors /3{6) = f3{6;n^ ,w) such that 

uniformly on compact subsets ofQ. 

This convergence also holds true up to the "nice" parts of the boundary; more- 
over, considering ratios of observables corresponding to different tn's, one can get 
rid of normalization issues. We work this out in Theorem 13.161 

A striking feature of our new observables is their direct relation to spin correla- 
tions. Let be a simply connected domain with m punctures, that is, m single 
faces /i, . . . , fm removed, and let w be the cover that branches around each of these 
punctures. Then, it turns out that F^{a,b), b € dD,^ , is (up to a fixed complex 



factor, see Proposition 3.6) equal to 

ZabKab[o-ifl) ■ ■ ■ Cr(/„i)], 

where Z^b and 'Eab stand for the partition function and the expectation for the Ising 
model with Dobrushin boundary conditions: "— " on the (ab) boundary arc and 
"+" on (ba). This, together with convergence results for the observables, gives the 
following corollary (with the notation "E+" referring to "4-" boundary conditions 
everywhere on dfl^): 

Corollary B. Let {Q^,a^,b^) approximate {n,a,b) as 6 0. Then 



E+[a{zf)...a{zf„)] s^o 



z?(0(zi),...,0(z„)), (1.4) 



where — i^q^q are explicit functions and 4> is a conformal map from f2 onto the 
upper half-plane C+ sending a to 00 and b to 0. 
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In Section |6] we give explicit formulae for d in C+, and hence, by conformal 
invariance, for all simply connected domains. For example, 

=cos [7rhmn(z,(a6))], 

where hmo(z, (ab)) stands for the harmonic measure of the arc (a6) in fl as viewed 
from z. These formulae for m = 1,2 were previously conjectured by means of Con- 
formal Field Theory, see |BG93| and earlier papers. To the best of our knowledge, 
the explicit formulae for m > 3 are new. 

Corollary [b] admits a number of generalizations. Let approximate a finitely 
connected domain Q with k inner boundary components 7i , . . . , 7fc (possibly macro- 
scopic). Then, for any m < k, one has 



C6(7l,---,7m), (1.5) 



where the functions '(^^^(TI; • ■ ■ j 7m) conformally invariant, the expectations Eq;,, 
E_|_ are taken for the Ising model with Dobrushin (respectively, ) boundary con- 
ditions on the outer boundary component and monochromatic on inner components 
7j, meaning that we constrain the spins to be the same along each component, but 
do not specify a priory whether it is plus of minus. In this case, tT(7|) denotes the 
(random) spin of the component 7^ . 

Further, closely following the route proposed by Hongler in [HonlOj . we prove a 



Pfafhan formula which generalizes (1.51 to the case of 2n boundary change operators 
(in other words, — / . . . / + /— " boundary conditions with 2n marked boundary 
points, see Section[5|. For m = 1, 2 this Pfaffian formula (along with the expressions 
for "d) was previously derived by means of Conformal Field Theory JBG93 , whereas 
we give it a rigorous proof for general m both in discrete, and, thanks to the 
convergence theorem, in continuous settings. 

Another application of our new observables |Izyll| is the proof of convergence 
of (multiple) Ising interfaces to SLE curves in multiply connected domains. In that 
context, a proper choice of the observable (i-e., the corresponding double cover 
w) guarantees its martingale property with respect to the growing interface. To 
prove that property, it is important to relate the values of to the partitions 
function of the model with relevant boundary conditions. In Section [Sj we show 
how to do it in the most general case, see Proposition [5^ The simplest example of 
an SLE process treated in this way (for which the use of a non-trivial double cover 
is essential) is a radial Ising interface converging to radial SLE3. 

For simplicity, in the present paper we work on the square grid, but all our 
proofs remain valid for the self-dual Ising model defined on isoradial graphs (e.g.. 



see [CS12j)- We refer the reader interested in a detailed presentation of the basic 
notions of discrete complex analysis on those graphs to the paper CSllJ and the 
reader interested in the history of the Ising model to the paper |CS12) and references 
therein. 

1.3. Organization of the paper. In Section [2j we fix the notations and conven- 
tions regarding discrete domains and the Ising model. In Section [3j we give the 
definition of the spinor observable and discuss its properties (in particular, discrete 
holomorphicity and boundary conditions), as well as the connections to spin corre- 
lations. We then define the continuous counterparts of the observables and briefly 
discuss their properties. Section |4] is devoted to the proof of main convergence 
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results for spinor observables: Theorem 3.13 (convergence in the bulk) and The- 
orem 3.16 (convergence on the boundary). We generalize our results to the case 
of multiple boundary change operators in Section [5j Finally, in Section [6] we give 
explicit formulae for the continuous observables in the punctured half-plane and 
for the scaling limits 'S^q appearing in Corollary [b| 
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2. Notation and conventions 

2.1. Graph notation. By (bounded) discrete domain (of mesh 6) fi* we mean a 
(bounded) connected subset of the square lattice (an example of a discrete 
domain is given on Fig. [l]) . More precisely, a discrete domain is specified by three 
sets V{n^) (vertices), J"(f2*) (faces) and £{n^) = Sini^l^) U £bd{n^) (interior edges 
and boundary half-edges, respectively), with the following requirements: 

• all four edges and vertices incident to any face / £ J^{n^) belong to £{fl^); 

• every vertex in V{il^) is incident to four edges or half-edges in 

• every vertex that is incident to at least one edge or half-edge e € £{V,^) 
belongs to V{n^); 

• at least one of two faces incident to any edge e G £in{^^) belongs to 7^(51"^). 

For an interior edge e S £in{^^) we denote by Zg its midpoint. For a boundary 
half-edge e G £hd{^^) we denote by Zg its endpoint which is not a vertex of i}^ . 
When no confusion arises we will identify an edge (or half-edge) e with a point Zg. 

By the boundary dfl^ of wc will mean the set of all its boundary half-edges 
£\yd{^^) or, if no confusion arises, the set of corresponding endpoints Zg. 

A double cover of a discrete domain 51* is a graph V,^ with a two-to-one local 
graph isomorphism w : (l^ ^ V,^ . Given a marked boundary half-edge a £ dfl^ , 
one can describe points z on a double cover by lattice paths 7 running from a to 2; 
in , modulo homotopy and modulo an appropriate subgroup of the fundamental 
group. If il^ is (fc-t- l)-connected, that is, has k holes, then there are 2*^ double 
covers, including the trivial one. Namely, to define a cover, for each hole one has to 
specify whether a loop surrounding this hole lifts to a loop in the double-cover, or 
to a path connecting points on different sheets. In the latter case we will say that 
vj branches around that hole. If z is a point on a double cover , we let z* £ fl^ 
be defined by n7(z*) — w{z) and z* ^ z. We will also use the obvious notation 
V(f2*), £{h^) etc. 
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2.2. Ising model notation. We will work with the low-temperature contour rep- 
resentation of the critical Ising model in fl^ (see |Pal07j ) . We call a subset S of 
edges and half-edges in (see Fig. [l] note that we admit inner half-edges in S) a 
generalized configuration or a generalized interfaces picture for this model, if 

• each vertex in V,^ is incident to 0, 2 or 4 edges and half-edges in S; 

• if an edge e — e'U e" consists of two halves e', e", then at most one of those 
three belongs to S. 

We will denote the set of all generalized configurations in fl^ by Confgcn(i^'^)- By 
the boundary dS of S* G Confgcn(f^'') we will mean the set of all half-cdgcs e £ S 
or corresponding points Zg, if no confusion arises. The partition function of the 
critical Ising model is given by 

Z{n^) = J2 :e = a;crit = \/2 - 1 (2.1) 

(the value x = Xcrit will be fixed throughout the paper). Here and below \S\ is the 
total number of edges and half-edges in S, and 

Conf(f^'^) := {S e Confgo„(l^'') : dS C dn^}. 



The formula (2.1) endows the set Conf(r2 ) of configurations, which corresponds 
to free boundary conditions in the spin representation, with a probability measure, 
the probability of a particular configuration S being a;'"^' /Z{n^). 

We will mostly work with subsets of Conf (fi"^), and restrictions of the probability 
measure to those subsets. Thus, we denote 

Conf+(f^'^) ~ {S £ Coni(n^) : dS = 0}, 

A A (2-2) 

where "mod 2" means that if some of ei, . . . , e„ appear several times in the sub- 
script (it will be useful for us to allow this), we keep in dS only those appearing an 
odd number of times. In the spin representation, the subset Conf+(ri'^) corresponds 
to locally monochromatic boundary conditions, that is, along each boundary com- 
ponent the spins are required to be the same (although they may be difl^erent on 
different components). If oi, . . . ,a2n & dfl^ , then Confai,...,a2„ (^'') corresponds to 
the configurations where the spins change from to "— " at the boundary points 
(half-edges) ai, . . . ,a2„. 

Remark 2.1. To simplify the notation, we will write Confej^^..._e„ (^^'') instead of 
Confro(ei),...,ro(e„)(^^'') whcu ci, . . . , e„ e ^l^. One should remember that we always 
consider Ising configurations or generalized interfaces pictures in the planar domain 
itself, even though we will define observables on double covers . 



3. Spinor holomorphic observables and their limits 

3.1. Discrete holomophic spinor observables. In this subsection we will con- 
struct spinor observables and prove their discrete holomorphicity. These observables 
should be considered as natural generalizations of fermionic observables introduced 
by Smirnov |Smi06[rCS12j to the multiply connected setup. A discrete domain f2*, 
its double cover tu : fi'' — >■ fi*, and a boundary half-edge a € dfl^ will be fixed 
throughout this subsection. In order to give a consistent definition for all discrete 
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Figure 1. An example of four-connected discrete domain 
and two generalized interfaces pictures Sa.b € Confa^z{fl^), each 
decomposed into a collection of loops and a simple lattice path 
"fA.B : a z, as required in Definition 3.1 For a double cover 
, branching, say, around two small holes but not around the 
big central one, one has 1{Sa) = 0, 1{Sb) = 2, and s(z,ja) = 
—s{z, jb) as, being lifted on , ja and jb end on different sheets. 



domains, we need the following (technical) notation. The half-edge a, oriented from 
an inner vertex to Za, can be thought of as a complex number. Then we set 

— g-f (arg(a) + i^) ^ (ia/\a\)-^ (3.1) 

for some fixed choice of the sign. Note that, if a is south-directed, then r/^ ~ ±1. 

Given a point z — z^ (i.e., the midpoint of an edge or the endpoint of a bound- 
ary half-edge e € E{n^)) and a configuration S € Confa^d^^), we introduce the 
complex phase of S with respect to z. First, we decompose 5* into a collection of 
loops and a path 7 running from 07(0) to 137(e) so that there are no transversal 
intersections or self- intersections (see Fig. [l]) . A loop in such a decomposition will 
be called non-trivial if it does not lift to a closed loop on the double cover (that 
is, if it lifts to a path between points on different sheets). We denote by 1{S) the 
number of non-trivial loops in S. We also introduce a sign 3(2,7) :— +1, if 7 lifts 
to a path from a to z on fJ"^, and 5(2,7) ■= ~1j if it lifts to a path from a to z*. 

Definition 3.1. The complex phase of a configuration S G Confj, ^(fi'^) with respect 
to a point z lying on a double cover ru : ^ V,^ is defined as 

W^{z, S) := e-4"(^) • (-l)'(^) • s(z, 7), (3.2) 

where w(7) denotes the winding (i.e., the increment of the argument of the tangent 
vector) along 7. Then, we define a spinor observable on the double cover il^ as 

F^{a, z) -.^ iija ■ W^iz,S)x\^\. (3.3) 

S6Conf„_j(n*) 

Remark 3.2. (i) Wt^{z, S) does not depend on the way one chooses the decomposi- 
tion of a given configuration S into loops and the path 7. The proof is elementary, 
and we leave it to the reader. Note that it is sufiicient to check that the second 
factor (— l)''^'^)s(2;, 7) is independent of a decomposition, as the rest is well known 
(e.g., see |HS10[ Lemma 7]). 
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(ii) By definition, F^{a,z*) = —F^(a,z), thus we call F^{a,z) a spinor. 

(iii) If 7Z7 is the trivial cover, then Definition 3.1 reproduces the original construction 
due to Smirnov (e.g., see |CS121 eq. (2.10)]). We denote this observable by Fo(a, z) 
and the corresponding complex phase by Wo{z, S). 



The most important "discrete" properties of the observable (|3.3|) are revealed 

Section 3] or 



in Theorem 3.3 below, which states its s-holomorphicity (see |SmilO 



|CS12( Definition 3.1]) and describes the boundary conditions. We introduce the 
following notation: given a vertex v £ V{fl^), we consider four nearby corners of 
faces incident to v, and identify them with the points Vk :~ v + e«'r(2's+i)/4 . 
fc = 0, 1, 2, 3. We denote sets of all cor ners of fl^ and its double cover O'^ by T(r2'') 
and T{il^), respectively. Similarly to (3.1), for a corner c — Vk £ T[^^) we set 
r;, := {i[c-v)/\c-v\)-^ := e-"(2fc+i)/8 

(again, the particular choice of square root signs is unimportant, so we fix it once 
forever for each of four possible orientations of c—v). We denote by 

Pr„^(F) Re(77^F) Vc = 1{F + vjF) 

the orthogonal projection of a complex number F Cz C onto the line rjc^. 

Theorem 3.3. For any corner c G T(f2'^) formed by edges or half-edges z' , z" € 
one has 

Pr„^(F^(a,z')) =Pr^,(i^^(a,z")). (3.4) 
Moreover, if b G dtl^ \ {a, a*} is a boundary half-edge, then F^{a,b) || rji,, i.e., 

Pr„,,(F^(a,6)) -0. (3.5) 

Remark 3.4. Since our observables satisfy F^{a,z*) = —F^{a,z), the identities 
(3.4) at two corners c, c* such that w{c) = izi{c*) are equivalent. The same is 
fulfilled for the boundary condition (3.5). 

Proof. Let v denotes the vertex incident to both z' and z". There exists a natural 
bijection 11 : Confa,z'i^^) — >■ Confa^z"{^^), provided by taking "xor" of a gen- 
eralized configuration S with two half-edges w(vz') and w{vz"). The well known 
proof of the theorem for the trivial cover (e.g., see |CS121 Proposition 2.5] or jHSlOi 
Lemma 45]) assures that, for any S G Confa,z'(^''): 

Pr,jWo{z',S)x\''\) = Pr,^(W^o(2",n(5))xin(^)l). 

Clearly, the same holds true with Wq replaced by Wro, unless 11 changes the number 
of non-trivial loops 1{S) or 5(2', 75) ^ s(z",7n(5)). However, it is easy to see that 
n always preserves the factor (—1)' • s: for instance, if there was a non-trivial loop 
in S that disappeared in 11(5'), then this loop has become a part of the path 7n(s)j 
leading to the simultaneous change of the sign s. 

To derive the boundary condition (3.5), it is sufficient to note that the winding 
of any curve 7 running from a to 6 is equal to (argfo — (arga-|-7r)) modulo 27r. □ 

The next proposition relates the boundary values of F^(a, •) to spin correlations 
in the Ising model. For a given double cover w : fl^ i}^ of a (/c + l)-connected 
domain fl^ , let us fix the enumeration of inner boundary components 71 , . . . , 7^ so 
that 

137 branches around each 0/71, ... , 7^ but not around 7m+ij ■ • • j 7fc- 
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For simplicity, below we assume that two marked boundary points a, b belong to 
the outer boundary of . We denote by Zab and Kab the partition function and the 
expectation in the Ising model with "— " boundary conditions on the counterclock- 
wise boundary arc {w{a)w{b)) C dQ^ , "+" on the complementary arc {w{b)TAj{a)), 
and monochromatic on all inner boundary components 71 , . . . , 7^ . Recall that these 
boundary conditions require the spin to be constant along each jj . We denote this 
(random) spin by cr(7j). 

Remark 3.5. If some hole 7^ is a single face, then we do not impose any boundary 
condition there and a{'^j) is just a spin assigned to this face. 

Proposition 3.6. If a,b E belong to the outer boundary and w{a) 7^ '^{b), 

then 

F^{a,b) = ±r]b ■ 'ZabKab[a-{ji)a{'^2) ■ • ■ ct(7„i)] (3.6) 



(the choice of sign is explained in Remark 3.1). In particular, Fo{a,b) — ±7/h • Zab- 
Moreover, 

F^{a,a) = irja ■ Z+E+[cr(7i)cr(72) . . .ct(7,„)]. (3.7) 



Remark 3.7. The sign ± in (3.6 1 depends on particular choices of 77^, rjb and the 
sheets of a, b on VL^. One way to fix it is as follows: let b be such that the counter- 
clockwise boundary arc {w{a)w{b)) C d^l^ lifts to {ab) C dV l^ (otherwise, consider 



F^{a,b*) ~ —F^{a,b)). Then, one can replace ±rib in (3.6) by —riaC 2*»''^ where 
Wab denotes the winding of the arc {-cu(a)zu(b)). 



Proof. The second identity is clear from the definition ( 3.3 ) , since each configuration 
S G ConfaJ^^) — Conf +(rj '^) contributes the same amount ir]a{—iy^^^x^^^ to both 
sides of (3.7|. To prove (3.61, note that each S € Confa^b{^^) contributes ±7]bX^^^ 



to both sides, thus we only need to check that all the signs are the same. Given a 
configuration S, decompose it into a path 7 : ■uj{a) Tu{b) and a collection of loops. 
A loop contributes to 1{S) (i.e., changes the sheet in Q^) if and only if it has an 
odd number of components 71, ... , inside. Hence, removing all those loops from 
the configuration results in the same sign change (— l)'^'^' for both sides. Removing 
the other loops (having an even number of those 7j's inside) does not change the 
signs. After removing all loops, moving 7 across any of 7j will change cr^jj) and the 
sheet on which the lifting of 7 ends (i.e., the sign 3(6,7)), again resulting in (— 1) 
factor at both sides. If finally 7 goes along the counterclockwise arc {■cu{a)nj{b)), 
then S contributes a;''^' to Za;,Eaf,[cr(7i)cr(72) . . . (7(7™)], while its contribution to 
the left-hand side is equal to i?7ae~5^^''''~'^-'a;l'^l = —riae~^'"''''x^^^ = ±?7f,a;l'^l, if 7 
lifted to the double cover ends at b, and riae^^'"'^''x^^^ , if it ends at 6*. □ 

3.2. Continuous spinors and convergence results. In this section we intro- 
duce continuous counterparts of the discrete holomorphic spinor observables, which 
we will later prove to be scaling limits thereof. For a moment, let us assume that $1 is 
a bounded finitely connected domain whose boundary components are single points 
71 = {wi}, . . . , 7s = {ws} and smooth curves 70, 7s+i, • . . , 7fc. Given a double cover 
ci7 : 51 — > £7 and a point a e 9r2\{wi, . . . ,Ws}, we denote by f^{a, ■) : ft C (or 
just for shortness) an analytic function which does not vanish identically and 
satisfies the following properties: 

(a°) /ro(a, z) = -fx^{a, z*) everywhere in fi; 
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(b°) /ro(a, •) is continuous up to dQ except, possibly, at the single-point bound- 
ary components and at a, and satisfies Riemann boundary conditions 

/ro(a, 2:)^/^r^^ e M, z e dil\ {a,wi, . . . ,Ws}, 

where denotes the outer normal to f2 at 2;; 
(c°) for each single-point boundary component {wj} the following is fulfilled: 
if w branches around Wj, then there exists a real constant Cj such that 

c ■ 

/ro(Q, z) = ^ ^ + 0(1) as Z~^Wj] 

otherwise is bounded near Wj, and thus has a removable singularity 
there; 

(d°) in a vicinity of the point a, one has 

/^(a, z) — ^' ° V 0{1) as z — > a 

z — a 

for a reaZ constant c". 
The properties (b°)-(d°) should be thought of as natural continuous analogues 
of those satisfied by F^{a,-) on the discrete level. Namely, (b°) corresponds to 
the boundary condition ( |3.5[ ); (c°) turns out to be the correct formulation of this 
boundary condition for microscopic holes; and (d°) states that has the simplest 
possible behaviour near a, which roughly resembles the fact that F^{a, ■) fails to 



satisfy (3.5) at one boundary edge only. 



Remark 3.8. Lemma 3.10 below shows that properties (a°)~(d°) define the function 
/ro uniquely up to multiplication by a real constant; moreover, c"^ ^ unless 
vanishes identically. Sometimes it is convenient to fix this constant so that 

c" ^ 1. (3.8) 

However, below we also work with non-smooth domains, for which c" is not well- 
defined; therefore, we prefer to keep the multiplicative normalization of unfixed. 

The boundary value problem (a°)-(d°) is not easy to analyse directly. However, 
the following trick relates it to a much simpler Dirichlet-type problem: given a 
spinor Z^, denote 



h^{z) :=Im J (/^(C))'dC, zen. 



Note that the function (/^(C))^ ~ (/ro(C*))^ is analytic in 51, so is locally 
well-defined and harmonic. 

Lemma 3.9. A holomorphic spinor satisfies the conditions (b°)~(d°) if and 
only if is a single-valued harmonic function satisfying the following properties: 

(bf^) ft-ro is continuous up to dn\{a,wi, . . .Ws}; moreover, = const and 
dnhru < on all macroscopic inner boundary components 7s+i, . . . ,7fc o-nd 
the outer boundary of Q, where i9„ stands for the outer normal derivative; 

(c°h) ^ro is bounded from above near single-point boundary components {wj}; 

fd'^) h^ is bounded from below near the point a. 

Proof. Let be a holomorphic spinor such that (b°)-(d°) are fulfilled. On the 
smooth boundary 9r2\ {a, wi, . . . , Wg}, one can reformulate the boundary condition 
(b°) as ifruiz))"^ -in^ > 0, which is equivalent to say that dinh^ = and (?„/i < 0, as 
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stated by (b^). Further, (c°) is equivalent to say that h{z) — Cj log \z ~ Wj\ + 0{l) 
as — > Wj, hence (c^) holds true. In particular, is single- valued in Q as it 
is single-valued near each of Wj and constant along each of macroscopic boundary 
components. Similarly, (d°) can be rewritten as h{z) = — (c'')^Re^3^ -I- 0(1) as 
z — >■ a, which is equivalent to (d^). 

Vice versa, if satisfies Dirichlet boundary conditions on smooth macro- 
scopic boundary components, then = [dyh^ + idxh^]^/"^ is continuous up to 
dfl \ {a, wi, . . . ,Ws}. Then, one can easily apply the similar arguments as above to 
deduce (b°)-(d°) from (b°)-(d°). □ 

Lemma 3.10. The holomorphic spinor f^{a, •) with the properties (a°)-((r ) above, 
if exists, is unique up to multiplication by a positive constant. Moreover, if (p : ^ 
n' is a conformal map, then, again up to multiplicative constants, 

f^{a,-) = {cl^'{z)Y f^,{<j>{a),cp{-)) (3.9) 

and h^{a, •) — /i^, (0(a), with vo' being the pushforward of m by (f>. 

Proof. Suppose that /i, /2 both satisfy (a°)-(d°). Then, one can compose a lin- 
ear combination / := C2/1 — 0^/2 with non-zero coefficients, such that / satisfies 
(a°)-(c°) and is bounded in a neighborhood of a. As above, define a harmonic func- 
tion h{z) :— Im / {f{C))^dC, and note that it is continous up to dil \ {wi, . . . , Ws}. 

As (b°) implies i9„/i < everywhere on macroscopic boundary components, the 
function h cannot attain its maximum there. Also, (c°) says that h is bounded from 
above in a neighborhood of each Wj. Then, the maximum principle gives h = const 
and / = everywhere in 17, i.e., fi and /2 are proportional to each other. 

For the second claim, it is sufficient to check the properties (a°)~(d°) for the 



right-hand side of (3.91 and apply uniqueness, which we leave to the reader. □ 



The existence of a non-trivial solution to the above boundary value problem 



will follow from Theorem 3.13 we also refer the reader to |HP12j for a purely 



analytic techniques developed for boundary problems of this kind. 



The conformal covariance property (3.9 1 immediately allows one to extend the 
definition of to non-smooth of unbounded domains: 

Definition 3.11. If is a finitely connected bounded domain, such that dVl con- 
sists of smooth arcs and single points, we define /^(a, •) as the unique, up to a 
multiplicative constant, non-zero solution to the boundary value problem (a° )~(d° ). 



Otherwise we define it by (3.9), where fi' is any smooth bounded domain. 

Further, we choose a harmonic function /i^(a,z) := / (/^(a, C))^'^C so that 
/i^(a, •) = Q on the boundary component of Vi containing the point a, thus h^ is 
defined up to a multiplicative constant as well. 

Remark 3.12. An equivalent definition of h^^{a, •) for non-smooth domains would be 



3.9 



to impose the conditions (b^)-(d^) given in Lemma 3.9 and the condition that f^ 
[dyh^^ -\- idxh.^Y^'^ is a spinor on f2. Indeed, the only condition in in Lemma 
that relies on smoothness of Sfi is dnh^^ < which can be reformulated in the 
following conformally invariant way: 

there is no b G Oil \ {a, wi, . . . , Ws} such that h.^{-) < /iro(&) near b. (3.10) 

For our convergence results, we assume that discrete domains approximate 
ri in the Caratheodory topology, see [Pom92j or [CSlli Section 3.2]. The reader 
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unfamiliar with that notion can think of the (stronger) Hausdorff convergence. To 
simplify notation, we also assume that i}^ has the same topology as f2. The first 
theorem says that discrete spinors defined in Section 3.1 (with respect to a fixed 
double cover w of the refining domains Q^) are uniformly close to their continuous 
counterparts in the bulk of fl. 

Theorem 3.13. Suppose that is a sequence of discrete domains of mesh size 
5 approximating (in the sense of Caratheodory) a continuous finitely connected 
domain ft, and that G dfl^ tends to some a £ d^l which is not a single-point 
boundary component. Then, there exists a sequence of normalizing factors (5(6) ~ 
l3{S;n^,a'^,w) such that 

uniformly on compact subsets ofQ. 

Proof. See Section |4] □ 

When extending this convergence to the boundary, we will impose additional 
regularity assumptions: 

Definition 3.14. We say that a sequence of discrete domains fl^ with marked 
boundary points b^ approximating a planar domain il with a marked boundary point 
b is regular at b, if 

• near b, the boundary dil locally coincide with a horizontal or vertical line; 

• there exist s,t > 0, such that, for any S, V{n^) contains a discrete rectangle 

R^{s,t) -.^ {S- {k + i{l + l)) : -s<k6< s, < 15 < t}, 

shifted and rotated so that b^ is the midpoint of its boundary side, and dil^ 
coincides with that side in the s-neighborhood of b^ . 

Remark 3.15. In fact, all our results can be directly extended to the case of a 
straight, but not necessarily vertical or horizontal boundary (cf. |CS12l Theo- 
rem 5.6]). Some additional technicalities are required to prove this result in the 
full generality: note that is not even continuous or bounded on the non-smooth 
boundary, so one is forced to work with ratios, as, e.g., in Theorem |3.16| below. 

Theorem 3.16. Let ■uJi,W2 be two double covers of a bounded finitely connected 
domain VL with two marked points a, b on the outer boundary component, and let 
57* converge to Vl in the Caratheodory sense, a^,b^ be boundary points converging 
to a, 6 as (5 — > 0, and this convergence is regular at a and b. Then, 

F^^{a\b')F^,{a\a') ^ /^,(a,6) 
where both /^i and f^^ '^^^ normalized by (3.8). 

Remark 3.17. Formally, above we should have used different notation 0^^,02 etc, to 
denote points lying on different double covers, with wi{a\) = 1172(02). We prefer 
to keep a simpler notation for the shortness. 

Proof. See Section |4] □ 

In the next corollary, let tu be a fixed double cover of fi, and 71, . . . ,7m be 
those inner components of for which w branches around 7j . Denote by 7^ the 
corresponding components of . 
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Corollary 3.18. Under the conditions of Theorem \3.1(^ as 5 0, one has 



where both and f^^ are normalized by (3.8). The limit '&abi'^ii ■ 
conjormal invariant of (f2; a, 6). 

Proof. Denote a{V) := cr(7f )cr(7^') 

Ea«b«[CT(r) 



(3.12) 

7m) is a 



(t(7^). Due to Proposition 3.6 

ZaabiEa-Sb^f (r)] Z+ 



one has 
F^{a\b^) Fo{a^,a^) 



E+[a(r)] Z+E+[a(r)] Z^.^. F^(a^a«) i^o(a^&')' 

Thus, ( |3.12 1 immediately follows from Theorem 3.16 The limit is conformally 



invariant due to the same conformal covariance property (3.9 1 of both and /q 



(observe also that the coefficients , eg for both and /o change by the same 
factor 10' (a) I -1/2 ^^^^ applying Q). □ 



Remark 3.19. (i) Corollary 5.10 gives a generalization of this result for the case of 
2?i+2 marked boundary points and "+/ — /.../ + /— " boundary conditions, 
(ii) Let a third point d be marked on the outer boundary of and the convergence 
of to n is regular at d as well. Then, 

Z,s,s ^ |fo(a^6^)| \fo{a,b)\ 
Z,s,s \Fo{a^,dS)\ |/o(a,d)|' 

and this limit is a conformal covariant of the multiply connected domain (fi; a, 6, d) 
(namely, it is multiplied by the factor \<l)' {b)\^/^\(j)' {d)\^^^^ when applying a confor- 
mal map (p). For simply connected J7's, this is given by |CS12i Corollary 5.7], and 
we give a proof for multiply connected domains in the end of Section [4] 

4. Proof of the main convergence theorems 



In this section we prove Theorems 3.13 and 3.16 following the scheme developed 
in |CS12j for simply connected domains. First of all, in order to transform the 
boundary conditions (3.5) to the Dirichlet ones, we consider a discrete integral 
Im /(F^(z))2d^TFurther, we observe that, under a proper normalization, 
the functions and have non-trivial subsequential limits on compact subsets 
of as (5 — )• 0. We then show that any such subsequential limit is a solution to 



the boundary value problem (a°)-(d°), and Lemma 3.10 guarantees that all those 



limits are the same, concluding the proof of Theorem 3.13[ Finally, we tr eat the 
behaviour of F^ near the boundary points b and a to prove Theorem 3.16 



For technical purposes, we extend our domain slightly: denote by dJ^{i}^) and 
dV{n^) the subsets of faces and vertices that are adjacent but do not belong to 
J^(r2*) and V{Q^), respectively. More precisely, dV{i}^) can be identified with the 
set £hd{^^) of boundary half-edges and should be formally considered as a set of 
pairs {{v;e) : e = (wintw),^ ^ V(f^'^),Wint e Vifl^)} (e.g., see [CSTTl Section 2.1]), 
and dJ-{il^) should be treated in the same way. We set 



V(fi*) := v{n^)udv{n^) 



and 



We work with the discrete spinor F{z) — F^{a;z) defined on a double cover 
of a discrete multiply connected domain $7'', which we do not include in the 
notation unless needed. Recall that it is s-holomorphic (i.e., satisfies (3.4)) and 
obeys the boundary conditions (3.5) at all boundary half-edges e, except for one 
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edge a on the boundary. We denote the corresponding vertex of dV{fl^) by Va- 
Observe also that is not identically zero, since the positivity of spin correlations 



and Proposition 3.6 yield [irja)''^ F^{a\ a) > 0. These are the only properties of 
F^(a; ■) that we will use in this section. 

Recall that, in the continuous case, it was proved to be useful to transform the 
boundary value problem (a°)-(d°) into a Dirichlet-type one by integrating The 
extension of this construction to the discrete setting is delicate, because the square 
of discrete analytic function need not be discrete analytic. Fortunately, the tools 



to treate this issue have already been developed in |Smil01 rCS12] . Proposition 4.1 
below summarizes these tools. Namely, properties (l)-(3) thereof show that one 
can define the discrete analog of Im J as a pair of functions and Ho , one 
of which is subharmonic and another superharmonic; properties (4)-(6) handle the 
boundary conditions, and properties (7), (8) show that these two functions cannot 
be too far from each other. This allows one to work with a pair H,, Ho as if it was 
a single harmonic function. Essentially, our analysis is based on a priori bounds for 
H derived from simple harmonic measure estimates combined with the uniqueness 
of solution to the boundary value problem (a°)-(d°). 

We define two functions H, and Ho on V{V,^) and J^{n^), respectively, by the 
following rule: if c e T{n^) and e G £{n^), v G V(n^), f € T{n^) are ah incident 
to c, then 

H,iv) - Hoif) V25 ■ |Pr,jF(Tn-i(e)))p. (4.1) 
Thanks to the square, this definition does not depend on the choice of the sheet. 



and thanks to the basic definition (3.4 1 of s-holomorphicity, it does not depend on 



the choice of e between the two edges (or boundary half-edges) adjacent to c. 

Proposition 4.1. The functions H, and Ho obey the following properties: 

(1) they are well-defined up to an additive constant; 

(2) ifee £(n^) is incident to fj' e and v,v' e V(f7'^'), then: 

H,iv) - H,{v') = \n^[{F{vo-\e)f{v - v')]- 
Hoif) - Hoif) - lm[{F{w-'{e)nf /')]; 

(3) AHoiv) < and AH,{f) > for all v G V{n^), f G J^i^l^), where A 
stands for the discrete Laplacian 

iAH)ix) ■.^^Y.iHix,)-Hix)y, 

(4) ifve dVin^) \ {va}, then H,{vint) - H,{v) > 0; 

(5) Ho = const —: Cj along each component ofdJ^{^l^) (then we fix an additive 
constant in the definition of H so that Ho = on the boundary component 
that contains a); 

(6) one can modify the discrete Laplacian at all vertices Wjnt G V{Vt^) incident 
to V £ dV{D,^) \ {va} and set values of H, on dV{il^) \ {va} to be equal 
to the corresponding Cj, so that ^ and will still hold (moreover, this 
modification do not destroy any estimates or convergence results for discrete 
harmonic functions defined on V{n^)); 

(7) if vi,V2,V3,V4 G V{^^) and fi, f2, fs, fi & J^{fl^) are adjacent to some 
inner face f G J^{il,^) and m := minHoifj), then, for all j = 1,2,3,4, 

iJ.(wj) - Ho{f) < const • (Hoif) - m) 
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with some universal constant; 
(^) */ /ii /25 /sj /4 G -^i^^) ^''^d Vi,V2,V:},V4^ G V{il^) are adjacent to some 
inner vertex v G V{n^) and M ~ max_ff,(wj), then, for all j — 1, 2, 3, 4, 



H,{v) - i/o(/i) < const • (Af - 
with some universal constant. 



Proof. All these properties are known in the simply connected case (e.g., see [CS121 
Section 3]). Since ([S]), ([T]) and ^ are local consequences of s-holomorphicity 
(3.4), they extend immediately to the multiply connected setup. Properties Q and 
(5j) follow from the boundary condition (3.51 using ([2|. The property ^ is also a 
local property of s-holoniorphic function satisfying the boundary condition (3.5), 
see |CS12|, Section 3.6] or |DHN1H Proof of Proposition 8] . So we only need to check 
([T]), i.e. that the summation of (4.1) along any loop gives zero. For homotopically 



trivial loops this follows from the local consistency of the definition (4.1 ) exactly as 
in the simply connected case, and extends to loops running around holes by ([5]). □ 

Remark 4.2. Note that we have an immediate corollary of properties ([l])-(|4]): if F 
is an s-holomorphic spinor satisfying the boundary condition (3.5) everywhere on 
dfl^ including the point a, then F = 0. Indeed, ([s]) implies that the corresponding 
function attains its maximal value on the boundary, and then = const due 
to the property iAi which now holds true everywhere on dQ^ . Moreover, similar 



arguments show that a solution to the discrete boundary value problem (3.4), (3.5) 



is unique up to a multiplicative constant (the proof mimics (CS12[ Remark 5.1 

In what follows, we denote by hm^a {v, 7) the discrete harmonic measure of a set 
7 C V{n^) in the discrete domain fi"^ viewed from a vertex v. Recall that hm^s {v, 7) 
is given by the probability of the event that the simple random walk starting at 
V hits 7 before dV{fl^) \ 7. We use the same notation hmf^s{f,'y) for the discrete 
harmonic measure of a set 7 C J^(il*) viewed from a face /. Essentially, we will 
use only the following elementary properties of hm*^ , which also are fulfilled for the 
discrete Laplacian modified near the boundary as it is mentioned in Proposition [IT] 
property ^ (see [CST2l Section 3.6] or [DHNlli Proof of Proposition 8]): 

• weak Beurling-type estimate: there exist absolute (i.e., independent of S, 
fi* and 7) constants C,p > such that the following estimate holds true: 

dist(2:;7) 



hm^,.(z,9r!'>\7) <C 



(4.2) 



where distQ* (z, K) means the smallest r > such that z and K are con- 
nected inside i}^ fl Br{z). 

uniform estimates for the exit probabilities in rectangles: if — R^{2s,t) 
is a discretization of the rectangle (—2s, 2s) x (0, t) with s > t, 70, 71 denote 
correspondingly the bottom and the top side of R^ , then 

Imz 



Ci • < hm^. (z, 71) < hm|. (z, OR' \ 70) < C2 



t 



(4.3) 



for all z £ R^{s,t), where Ci, C2 > are some absolute constants. 



Note that (4.2) and (4.3) hold true for arbitrary isoradial graphs (see |CSllj ). 

Now w e proceed to the definition of normalizing factors (3(6) used in Theo- 
rem 3.13 Let a sequence of discrete domains i}^ approximate a finitely connected 
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Figure 2. (A) The proof of the estimate < const • /Jj,. in 
Lemma [44} If —Ho is big at /, it is also big at vertices near some 
path ^ : f Wj. Since the harmonic measure of 7 seen from v 
is uniformly bounded from below, — (w') has to be big as well. 
(B) Similar considerations near a give the estimate < const ■ (3^. 



planar domain il, whose boundary dfl consists of single-point inner components 
71 = {wi},...,7s = {ws} and continua 7s+i, . . . 7fc,7fc+i, where 7^+1 denotes the 
outer boundary oJVL. Let > be chosen sufficiently small so that, for any r < r*, 

has the same topological structure as fi, where -Br.si(a) denotes the proper con- 
nected component of Br{a)r\^l. Further, let 5* = <5*(r) > be chosen small enough 
so that, for any <5 < (S*, one has — a\ < \r and 7^ C Bi^{wj) for all j = 1, .., s. 

Let ^ F^{a\ • ) be the spinor observable (fs^sl in Vt^ and Hl,Hi be the corre- 



4.1 



sponding discrete integrals Im J {F^{z))'^d^z defined by (4.1 1 and Proposition 
We introduce normalizing factors > by 

(/3,*)-i :=maxf2.(,,)|i?*|, (4.4) 

where, similarly to the continuous setup, ri^{r) := 57* \ [Uj^]^i3*(u'j) U i?*Q(a)] and 
Bf.{wj),Bf.^{a) stand for discrete r-neighborhoods of wj and a in 57*. 



Remark 4.3. (i) Note that the right-hand side of (4.4 1 does not vanish. Indeed, if 
it did, F^{a^ , •) would also vanish identically. 

(ii) Since the limiting function has singularities at a and Wj , it is natural to cut 
them off when taking a maximum to get a correct normalizing factor. 

This choice of guarantees that |/3*if*| < 1 in 57*(r). Our next goal is to show 
that, for a fixed r, the functions fi^H^ are uniformly bounded on any compact 
subset of n. This is achieved by the following lemma: 

Lemma 4.4. For any r < \r^, the ratio (32rl Pr bounded uniformly in S < S^{r). 
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Proof. Denote n\r) := n^(r) \ B^^ ^{a) and 

By definition, n^{2r) C h^{r) C rJ-^Cr), hence (3^ < M < ■ We first show that 
Pirl is uniformly bounded. Recall that is subharmonic, is superharmonic, 
Hi{v) > H^{f) for adjacent v and /, and = on the boundary component 
containing a. Therefore, 

(Pl)-^ = max |maxg^,(^) ; - mingj-^,^^^ ijf | . (4.5) 

Note that in fact 



max^ 



^an^r) H, < max9o^(2-r) H, . (4.6) 

Indeed, the function is subharmonic in B2j,i'Wj)r\il^ , j = 1, . . . , s, and it cannot 
attain its maximum on the boundary component 7^ C B2r{wj) because of property 



(4) of Proposition 4.1 (and, if 7* consists of one face only, then is subharmonic 
everywhere in B2r{wj)). 



Thus, either (Z?*)^^ < iPir)^^ ^^'^ then there is nothing to prove, or 
(^^ ) - 1 = -i/f (/) for some / G dB'^ {w, ) . 

Since is superharmonic, in this case there exists a path of consecutive neighbors 

7 - {/ = /i ^ /2 ^ . . . } such that -{Mr' = H^ifi) > Htih) >■■■■ This path 
can only end up at 7^, where the superharmonicity of fails (see Fig.[2|\). Denote 
by 7' the set of vertices adjacent to 7. Then, the property (8) in Proposition 4 
ensures that 

h! + 0tr' < const • ((/3t)"^ - Hi) everywhere on 7', 

where we have used the estimate max^« (wj)^^ < (/?2r)~^ which was explained 
after (4.6). Hence, for some absolute constant pi > 0, 

Hi < -pii^)-^ + everywhere on 7'. 

Further, there exists a constant P2 > independent of r and S such that, for 
any vertex v e dB^^iiJUj), one has hm(i;,7') := liin^5(2r)us| («)j)(^'7') > P2 (see 
Fig.[2]A_). By subharmonicity of H^, this implies 

Ht{v) < {-pMr'+ iPlr') ■ hm(t;,7') + {Ply^ ■ (l-hm(i;, 7')) 

Since Hl{v) > -(/3i)~\ we infer that P^JPf. < 2(pip2)"^ 

It remains to prove that (3^/f3f. is uniformly bounded, which is done by exactly 
the same argument with the roles of and interchanged, inequalities reversed, 
and a playing the role of Wj . Note that the uniform lower bound for the harmonic 
measure hm(/, 7') still holds (this time, 7' will be the set of faces adjacent to a 
vertex path 7 which terminates at a, see Fig.[2|3), provided that / is chosen, say, to 
be the closest face to some fixed z e f^\-B2r,n(a) (e.g., see [CSlli Lemma 3.14]). □ 

Now we are able to claim precompactncss of the families {H^}, {F^} as J — ?> 0. 
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Lemma 4.5. Fix some sufficiently small r > 0. Then, there exists a subse- 
quence (5 = — )■ such that the functions j3^H^ converge to a harmonic function 
if : — > M uniformly on compact subsets of O. Moreover, the functions I3^{F^)'^ 
converge to := dyH + idxH uniformly on compact subsets offl. 



Proof. Definition of (3^ and Lemma 4.4 guarantee that P^H are uniformiy bounded 



in fi*(2 ''r) for any fc > 0, and hence on all compact subsets of ft. Due to |CS121 



Theorem 3.12], the functions pfH^ and the spinors ^J~^F^ are thus equicontinuous 
on compact subsets of ^l^ , and the Arzela-Ascoli theorem implies their subsequential 
convergence to a continuous function H and a spinor F, respectively. Morera's theo- 
rem, together with the discrete holomorphicity of F^ ^ implies that F is holomorphic. 
Since increments of j3f.H^ are given by discrete integrals Im [[3^.{F^ {z)Y z], one 
has H — Im J so ii is a harmonic function defined in Vl. □ 



The next step is to show that all these subsequential limits solve the correct 
boundary value pro blem . It is convenient to work with the function H . Recall that, 
(properties (|5j) and (|6|), = Cj on each of macroscopic 
By definition (4.4) of the normalizing factors 



4.1 



due to Proposition 
boundary components j[ 

f3f, we have < 1. Thus, taking a subsequence once more, we can assume 

that 



s+l: ■ 



c-i such that |c 



for some constants 



J — ?> Cj as S 



^0, j = s + l, 
< 1. 



,k + l, 



(4.7) 



Lemma 4.6. For any subsequential limit H from Lemma 4-5 satisfying {4-'"^h the 
conditions (b'^)~(d'^) from Lemma 



3.9 



hold true. Moreover, supfjj^j \H\ = 1. 



Remark 4.7. Below we use the equivalent reformulation (3.101 of the boundary 
condition dnH < which does not rely on the smoothness of dft. 

Proof. Property (b'^). Our first goal is to prove that H satisfies the Dirichlet bound- 
ary conditions on all macroscopic boundary components jj, j = s+1, . . . , k-\-l. Fix 
some small r' > 0, and recall that IP^H'^l < Cir') everywhere in U,^{r'), where C(r') 
does not depend on 6. Due to superharmonicity of H^, for any face / G 51* (r'), one 
has 

PtHiif) > ■ hmf,.(,,)(/,7^) - C7(/)(l-hn4(,,)(/,7|)) ■ (4.8) 

Similarly, subharmonicity of implies that, for any vertex v e f2'^ (r'), one has 

< Ptcl ■ hmf,.(,,)(«,7|) + C(r')(l-hmf,.(,,)(z;,7f)) . (4.9) 

Now let / and v approximate some point z G ^{r') as (5 — > 0. Since both discrete 
harmonic measures converge to the continuous one, and H^^{f) < H^{v) for incident 
and ( |4.9P yield 

Cj hmn(,,,)(z,7j)| < C(r')(l- 



/ and V, (4.< 



\Hiz) 



-hm 



In particular, H{z) — cj as z tends to 7j. The boundary condition (|4| from 



Proposition |4. 1| survives in the limit and gives (3.10) due to |CS12l Remark 6.3]. 

Properties fc°J and (d'^). By subharmonicity of H^, the inequality f^^H^ < 1 
extends from dB^{wj) C dfl^{r) to the whole discrete disc Bf.{wj), thus 7i < 1 in 
a vicinity of Wj. Similarly, superharmonicity of and the condition = on 
the boundary component containing a imply H > —1 near a. 
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Normalization sup^^^) \H\ = 1. Recall that, similarly to (4.5), \(3f.H^\ attains its 
maximum max^^s^^) \/3f.H^\ = 1 on the boundary 

dn'ir) C a<n(a) U dB^w,) U • • • U dB'.iws) U 7^+1 U • • • U 7^+1. 
Note that (|4.7p and the property (b°J yield 



Moreover, as P^H^ uniformly converge to H on all compact subsets of ft, one has 

ma,-XQB^{wi)u---udB^(w^) ll^^H^l — > raa.xgB,^(wi)u---uaBr(w^) \H\. 
Finally, the convergence of fi^H^ to H inside of also imply 



since the estimates (4.8), (4. 9) guarantee that \I3^H \ are uniformly small near the 



boundary component of f2 containing a. □ 
Proof of Theorem 13. 13[ We set I3{5) := for a smaU fixed r > 0. By 



Lemma 4.5 the functions {(3{S))'^H^ and I3{S)F^ have subsequential limits and F, 
respectively. Lemmas [4.6[|3.9| and |3.10| guarantee that all these possible limits H are 
the same. Mor eove r, this unique limit is nontrivial due to the normalization condi- 
tion in Lemma 4.6 and thus coincides with /i^, normalized so that sup^^^) = 1. 
Since H{z) = {F{C))^dC, we conclude that F = f^. □ 

The next lemma shows that the convergence of f3{S)F^ holds true at the boundary 
point b (in fact, it follows from our proof that this convergence is uniform up to 
straight parts of the boundary d^l). Note that a similar result has been obtained in 
|CS12[ Theorem 5.6] under milder assumptions. For convenience of the reader, we 
give a shorter proof here using our (stronger) regularity assumptions for near b. 

Lemma 4.8. Suppose that, under the conditions of Theorem \3 . 1 3[ we are given a 
sequence of marked points b^ £ dQ^ , b^ — > b, and the convergence Q'' — )> £7 is regular 
at b. Then, l3{S)\F\b^)\ |/S(6)| as S 0. 

Proof. Below we assume that il and fl^ are shifted and rotated so that b^ — b ^ 0, 
il^ contains a discrete rectangle R^{s,t) for some fixed s,t > 0, and dil,^ locally 
coincides with the boundary of the discrete upper half-plane (see Defini- 
tion 3.14). We also assume that the additive normalization of the functions is 
chosen so that they vanish on the macroscopic boundary component containing b. 

Recall that, on compact subsets of fi, f3{6)F^ and {[3{S))^H^ converge to the 
properly norm a lized functions and h^, respectively. Moreover, the uniform 
estimates (4.8), (4. 9) guarantee that the convergence {/3{6))'^H^ — )• /i^, remains true 



up to dil, \ {a,wi, . . . ,Ws}- In particular, the functions {^{Syf'H^ converge to /i^ 
uniformly in the fixed rectangle R{s,t) around b. 
Let 

M := (/^(O))' - dyh^iO) 

and 



ntif) ■■=lmf, for faces /e 

s 

V2 



Hiiv) := Imv + -75 , for vertices v e C;^ 
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(note that Ti.^ can be defined by (4.1 1 starting with th e con stant s-holomorphic 
function = 1 which satisfies the boundary conditions (3.5 1 on 3C']_). Then, for 
any e > 0, one can find a small d > such that 

\{P{5)fH^ - ^M^\ < ed everywhere in := R^{2d,d) 

for all sufficiently small S^s. Since both Ho and Hi are discrete harmonic and 
satisfy the same Dirichlet boundary conditions as near , the sub- and super- 
harmonicity of and implies that 

Imrn^, - f^ni < ed ■ hm^. ( . , dR' \ dC%) ^-e'^y-here m R . (4.10) 

Let :— be the inner vertex of the boundary half-edge h (see notation on 
Fig. SK.). Estimating the discrete harmonic measure hm^s {vb, dR^ \ dC\_) by (O), 
one obtains 

V2d ■ |/3((5)i^*(6)cosf p = H^,ivb) < i^nt {v b) + const ■ e6. 
Since nHvb) = f + 73 = (\/2cos2|)5, this gives 

\l3{S)F\b)\^ < ^ + const -e. 

Note that this bound also holds true for the value \I3{S)F^ {b')\'^ , where b' :— S 
denotes the neighboring boundary edge and :— (1 -I- ^i)d is the corresponding 
inner vertex (see notation on Fig. [sj^). Now let b" := + ^i)S be the midpoint 
of the edge (wf,w^), and fb '■= (| -I- i)S be an inner face incident to both Vb and 
Using (4.101 and estimating hmj^s{fb,dR^ \ dC^_^) by (4.3 1, one obtains 

Re{l3{6)F\b")f ■ 5 = H^oih) > - const • e6. 
Note that for any complex number ^ G C one has 

Reie) < [2cos2f]-i • (|Pr,-ws(OP + |Pr,>,/B(u;)n. 
Therefore, we arrive at the inequalities 
H - const • e < Re{/3(S)F\b")f 

< [2cos2f]-i • (|Pr,.„/B(/3(5)F^(6"))P + |Pr,„/s(/3(5)F*(&"))n 
= [2cos2f]-i . (|Pr,-„/.(/3(J)F*(6))|2 + |Pr,,,/«(/3(<5)F*(6'))n 
= i(|/3(5)F^(6)|2 -I- \fi{6)F\b')\^) <^l + const • e. 

Since e can be chosen arbitrary small, this yields lim5_j.o \ fi{5)F^ {b)\^ = ^. □ 

Finally, we work out the relation of the values of discrete observables F^ at the 
points to the growth rate of their limit near a. 

Lemma 4.9. Under the conditions of Theorem \ 3.13\ assume also that the con- 
vergence fl^ ^ n is regular at a. For j — 1,2, let Wj be double covers of fl^ , 



7S 
TO 

f^. near a (see (dP ) in the definition of f^- given in Sect. 3.2). Then. 



l3j{5)F^, — /roj, and c°j denote the corresponding coefficients in the expansions of 



MS)Fi^ cj 

/32(<5)F4(a^) ci ^ ■ > 
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rf< 1 




R{2d\ d) 




oh 


Vb' 
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{A) 



b' 



R{s,t) 




Figure 3. (A) Notation near used in the pro of of Lemma 
(B) Notation near used in the proof of Lemma 



4.9 



Proof. Below we assume that and fl are shifted and rotated so that a-^ = a = 0, 
51* contains a discrete rectangle R^{s,t), and dfl^ locally coincides with dC^. We 
also write for F^^ ^ and /i^2 for /roi,2- Denote 

Recall that both F^ia) are positive multiples of some fixed complex number ii^a 



(see Proposition 3.6 1, thus > 0. Taking a subsequence, we may assume that 
^ k < +00 as (5 (if fc = +0O, then swap F^ and F^ and consider the 
inverse ratio (K^)^^). Note that 



f3US)FUz) ^ /i,2(z) 

(5^0 



0(1) inside i?(s,t). 



Therefore, it is sufficient to prove that the function 

F\z) ■ p^iS) Fiiz) ~ (5,{S) Ff{z), 

defined in R^{s,t), converges to a limit which remains bounded near a = 0, since 
this will immediately give — = for any subsequential limit. 

Note that, being a real linear combination of Ff and F2, the function F^ is s- 
holoniorphic in the discrete rectangle R^{s,t) and satisfies the boundary condition 
(3.5) on its bottom side, including the point a*, where F^{a^) = 0. Thus, (4.1) 
allows one to define a discrete integral Im {F^ {z))'^d^ z inside R^{s,t) so 

that = everywhere on the bottom side of R^{s, t). 

We claim that both functions f3j{6)Fj, and hence F^ and , are uniformly 
bounded on the top, left and right sides of the smaller rectangle R^{^s,d) (see 
Fig. [3^ for the notation), where d s. On the top side, this follows from the 
uniform convergence of (5j(S)Fj to a continuous limit. In order to prove a uniform 
bound on the left and the right sides, note that, for u e ± ^s + R{2d, 2d), the second 
terms in (4.8) and (4.9) can be uniformly estimated using (4.3) in the following way: 

1 — huifisf^r') ("j 7o) — const (s, t, r') • Imit. 

Therefore, {/3j(6))^H^{u) = O(Imu) in nei ghborhoods of the left and the right sides 
of R^{^s,d). Hence, l3j{S)Ff = 0(1) on those sides due to j(7S12l Theorem 3.12]. 

Thus, is uniformly bounded on the top, left and right sides of the rectangle 
R^{^s,d) and vanishes on its bottom side. Using super-/sub-harmomcity of 
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on faces/ vertices, and uniform estimates (4.3), one easily deduces from here that 
H^{u) — O(Imw) everywhere in R^{^s,d). Applying |CS121 Theorem 3.12] once 



0. 



□ 



again, we conclude that ~ 0(1) in R^{^s. d). Hence, — 
Proof of Theorem 13.161 One has 

F^M')F^Ab') f3i{S)F^M) ■ /32{S)F^,{b^) ^ (c^)-i/g, («, ' 
where we have applied Lemmas |4.9| and |4.8| to the first and the second factors. □ 



Proof of Remark 3^9](ii). Let r > be chosen small enough. As above, it is 
sufficient to prove (3.131 for any subsequence S — 6k such that /3{S)Fq converge 



to some nontrivial continuous limit. But this immediately follows from Lemma |4.8| 
applied to both boundary points b and — > d, since the properly normalized 
observables /3((5)i^Q (a'^, • ) converge to f^{a, ■ ) at both b and d. □ 



5. Multiple boundary change operators 

In this section, we follow |HonlO| to extend the definition of the spinor observ- 
ables to the case of multiple marked points on the boundary. We show that these 
observables (we call them multi-source ones) are still s-holomorphic. Moreover, by 
analysing boundary value problems they solve, we prove recurrence relations that 
eventually allow one to express all these observables in terms of the basic ones 
introduced in Section [3l 

The reason to introduce the multi-source observables is revealed in Propositions 
|5.4| and |5.6[ where we establish their relation to spin correlations and partition func- 
tions, respectively. The latter is especially important in view of two applications. 
First, it allows one to prove the discrete martingale property of those observables 
with respect to interfaces growing in multiply connected domains, leading to a de- 
scription of scaling limits thereof |Izyll| . Second, in the case of 2n microscopic 
holes carrying one boundary change operator each, the Kramers- Wannier duality 
relates the corresponding partition function to the 2n-points spin-spin correlations 
in the critical Ising model with free boundary conditions. This can be used to 
prove the conformal covariance of their scaling limits, which would complement the 
results of (CHI12J and this paper. 

Consider a domain with 2n-|-l marked points a :~ ao, ai, . . . , a2n G dfl^ 
and an inner edge z. Each configuration S € Coiiia(,^ai,...,a2„,z can be decomposed 
into a collection of (mutually disjoint, non-self-intersecting) loops and n + l curves 
connecting a^'s and z in some manner. In order to define the complex phase of S, 
we draw n artificial arcs vi, . . . , Vn connecting w{ai) to 07(02), • • ■ , ^(a2n-i) to 
'Cf7(a2n), respectively, and fix the way how they lift to the double cover 57. Adding 
these arcs to a configuration promotes it to a collection of loops and a single curve 
7 running from 'cc7(ao) to w{z). As we admit intersections of the artificial arcs with 
curves constituting 5, this time the loops can be self-intersecting, see Figure [4] 

Definition 5.1. We define the complex phase W^{z, S) to be e^s^^'''' (— l)'^'^'s(2:, 7) 
times (— 1)"'^'^'^^ where I{S) is the number of loops in SUviU ■ ■ ■Ui'n that have zero 
winding modulo An, and other factors are as in Definition \3.1\ Further, we define 
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Figure 4. (A) A doubly connected domain il'^ and a configura- 
tion S £ ConfaQ,ai,...,a6;z{^^)- Adding artificial arcs (dashed lines), 
we promote S to a collection of loops and a simple path 7 running 
from w^ao) to w{z). The loop containing 05 and ag has winding 
0, thus contributing to I{S). (B) A configuration contributing to 
/(ao, ai, a2, a3, 04, 05, flg; 03). It also contributes the same value 
(up to a complex factor e~i^^''^-°'^^°'^^) to /(oq, ai, 02, 05, 05; 04). 



the multi-source observable by the formulae 
F^{ao,ai, . . . ,a2n;z) iT]a ■ 

SeConfa 



E 



\s\ 



(5.1) 



Remark 5.2. In fact, the only data we use concerning each of the artificial arcs is 
its winding w(vs) modulo An and the way how it lifts to fl. The reader may check 
that altering the choice of {I's} can only result in a sign change of the observable. 



The following straightforward generalization of Theorem |3.3| holds true: 

Proposition 5.3. The observables {5.1) are s-holomorphic spinors satisfying the 
boundary condition (3.5) everywhere on d(l^ \ {ao, . 



Proof. The proof is essentially the same as the one of Theorem |3.3| The only differ- 
ence is that now the bijection 11 between the sets of configurations Confa,ai,....a2„,z' 
and Confa.ai,...,a2„.z" can create or destroy a loop which is self-intersecting. If such 
a loop has winding 27r modulo 47r, then there is no difference with the case of a 
simple loop. If it has winding modulo 47r, then its contribution to e~i^^'^^ after 
it becomes a part of 7 is minus that of a simple loop, which is compensated by the 
simultaneous change of I{S) by one. The boundary conditions follow in the same 
way as before. □ 

We now turn to the relation of the multi-source observables with spin correlations 
and partition functions. Let 2n + 2 distinct boundary points a = Oq, Oi, . . . , a2n, 
0-211+1 = 6 be chosen on the outer boundary of VL^ ^ and let to be a double cover that 
branches around boundary components 71 , . . . , 7„i and does not branch around the 
others. Write also :— cr(7i) . . . (7(7^), with aiV) = 1 if m = 0. 



Proposition 5.4. We have 



(5.2) 



24 



D. CHELKAK AND K. IZYUROV 



where the subscripts refer to the model with boundary change operators at all the 
marked points a, ai, . . . , 02^, b. 



Proof. The proof follows from definition of (a, A; •) similarly to (3.61, so we leave 



it to the reader. □ 



Remark 5.5. Let us describe one way to fix the sign in (5.2). Suppose that, starting 
at a and tracing the outer boundary of 17"^ in the counterclockwise direction, we 
find the marked points a, ai, . . . , 02,1, b in that order. Let the artificial arcs i^s run 
outside the domain, following the boundary arcs (02^-1, a2s). Then, assuming that 
the arcs (CT(a2s-i)cp(a2s)) car ry "+" boundary conditions, one can replace zLrji, 



in the right-hand side of (5.2 1 by —rjaC 2™'^''^ where Wab is the winding of the 



counterclockwise boundary arc {w{a)n7{b)), cf. Remark 3.7 



For the next proposition, we allow the marked points a — qq, . . . , a2n+i = & to be 
anywhere on the boundary. Denote by wz the double cover that branches around 
all boundary components of carrying an odd number of marked points and does 
not branch around the others. 

Proposition 5.6. Let be a discrete domain. One has 

(a, A; b) = ±776 • Za,A.b- 

Proof. Any configuration S € Confa.A,b contributes the same value irjbX^^^ to both 
sides of the equation, so it is sufficient to check that W{b, S) does not depend on S. 
It is convenient to add one more artificial arc connecting b to a, and to think that 
the arcs i^O: ■ ■ ■ :>^n are actually drawn on 57* as simple, mutually non-intersecting 
curves, cf. Remark [5^ 

By definition, an inner boundary component has an odd number of artificial arcs 
Vs issuing therefrom if and only if wz branches around that component. Hence, 
cutting along, say, the right-hand side of each Vg, we get a sheet of , and can 
be obtained from two copies of that sheet by gluing along the cuts, in such a way 
that crossing any cut would change the sheet. 

We now want to show that W{b,S) = e-5'^('>')(-I)'(-5)+^('5) 5(6,7) does not de- 
pend on S. Adding i/q to 5 U i^i U • • • U i^n, we end up with a collection of loops; 
denote 70 := 7 U i^q. Observe that a loop has winding modulo 27r (and hence 
contributes to I{S)) if and only if it has an odd number of self-intersections, and it 
contributes to l{s) if and only if it intersects an odd number of cross-cuts. So, we 
can write 1{S) + I{S) — X^ioops 70^70 ''(7") mod 2, where r(7a) is the number of 
intersections of \ {i^i, ■ • ■ ,J^n} with other loops (self-intersections contribute to 
both 1{S) and I{S)). Those intersections can only occur between S and artificial 
arcs, and each intersection contributes at most once to the sum. 

Further, observe that e~i""'^'^^ = ei""'^''°^e~i""^'^°K Just as above, -e-J'^^T") 
counts the number of self-intersections of 70, and 3(6,7) describes the number of 
intersections of 7 with the other cross-cuts. Therefore, 

W{S,b) = -e5™('''')(-I)'^*'^°'+^i°°p=^c.?^7o''('>'°). 

Since the exponent is the total number of intersections between the loops (not 
counting self-intersections), it is always even. □ 



In accordance with our convention "mod 2" in (2.2), Definition 5.1 also gives the 



values of at marked boundary points. Denote A :— {ai, . . . , a2„}. 
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Proposition 5.7. The identity 



2n 



F^{aQ,A;z) = ^ 



fc=0 



F^{ao,A;ak) 



(5.3) 



is fulfilled for any z £ . 



Proof. Both sides of (5.3) are discrete s-holomorphi c sp inors (defined on the same 
double cover f2*) satisfying the boundary condition (3.5 1 everywhere on dfl^ except 
the marked points oq, . . . , a2n- Moreover, for any fc = 0, . . . , 2n, there is only one 
term in the sum (5.31 that fails to satisfy ( |3.5| at a^. However, its value at coin- 
cides with the left-hand side value F^(aQ, A; a^). Hence, these two s-holomorphic 
spinors are equal to eac h oth er due to Remark [42} since their difference satisfies 
the boundary condition (3.51 everywhere on dH,^ . □ 



Note that configurations contributing to F^(aQ, A; a^) actually have 2n boundary 
points instead of 2n+2. Thus, the right-hand side of (5.3 1 can be expressed in terms 



of the similar observables with smaller number of marked points, and, recursively, 
in terms of the basic observables i^ro(afc;-). In order to do this in a convenient 
way, we need an additional notation. Recall that, for each a^, we fix the complex 
number rjk rja,. according to (|3.1|). Then, we define the real antisymmetric 



(5.4) 



(2n+l) X (2n + l) matrix = [(Gro)j.fc] by setting, for < j < fc < 2? 



(G^) 



j-k 



F^{af,ak) _ F^{ak;aj) 
iFr^{ak;ak) iF^{aj;aj) 

where we have used that {irikYe'™'^'^''^^ = 77^ for any curve ^kj '■ ^{o-k) ~^ nj{aj) 
running in fl^ . Further, let G^[ki, . . . ,ks] denote the sub-matrix of obtained 
by removing rows and columns with indices /ci, . . . , fcg. 



Proposition 5.8. We have 

2n 

F^(ao,A;z) = ±^(-l)'=PfG^[fc] •f^^(afc;z). 



(5.5) 



fc=0 



with the sign depending on the choices made for rjk and Vg ■ In particular, the sign 
is "plus" with conventions described in Remark \ 5. 9\ below. 



Remark 5.9. (i) The sign in the left-hand side of (5.5) depends on the choice of 



artificial arcs i^s, while for the sum in the right-hand side it depends on the sheets 
of ttk and the signs of r/fe, fc = 1, . . . , 2n. We will assume that each lifts to a path 
from a2s-i to a2s on 17*, and that the signs of rjk are chosen so that 



exp 



-w(z/s : a2s-i a2s) 



(5.6) 



(ii) One could also write (5.5 ) as the Pfafhan of a (2n-f 2) x (2n+2) matrix, obtained 



from Gro by adding the last column with entries F^ [ok ; z) and a corresponding row. 

Proof. We prove the claim by induction in n, starting with the trivial case n — 
and using (15. 31). For fc = 1, . . . , 2n, let k' := k + 1, if k is odd, and k' := k — 1, if 



k is even. By definition, any configuration contributing to Fj^{aQ, A; a^) contains a 
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curve running from to a^' , appended with an artificial arc connecting Uk' to ■ 
Removing tliat arc, and taking into account ( |5.6[ ), we get the foUowing identity: 

F^iao,A;ak) - {-!)'' ^F^{ao, A[k, k']; ak'), 

where A[k, k'] := A\ {a^, a^'}- Thus, using the induction hypothesis and observing 
that i^k^FlukTak) does not depend on fc, we get (for fc = 1, . . . , 2n) 



F^{ak\ak) irjk' 



J2 {-iypfG^[j,k,k'] 



0<j^k,k'<2n 



,1 F^{aj;ak') 
F^{ak]ak) 



, , F^{aj;ak' 



0<j#fe,fc'<2n 

Due to the standard recursive formula for Pfaffians applied to the matrix G^[k]^ 
this can be written as 



F^{ak;ak) 



(-l)'=PfG^[fc], /j = l,...,2n. 



(5.7) 



Similarly, 



irio Fro(a2, A[l,2];afe) 



i^tj,(ao,A;ao) = — Fro(a2, A[l,2];ai) = — E 

k—o 



F^{ak; ai). 



Applying (5.7) and using yyp F^{aQ;aQ) — rj^ F^{ai;ai), one arrives at 
F^{ao,A;aQ) _ ^ ^" 

fe=3 



F^{ao;ao) 



Plugging this and (5.7) into (5.3), we obtain (5.5). 



□ 



Corollary 5.10. With the notation as above and conventions of Remark 5.5 let 
57* approximate il as S 0, regularly at all marked points. Then, 

— 1 — ^ , (.5-oj 



E+[a(r)] 



Pf [Ca jak ] 0<j<A;<2n+l 



where the conformal invariants are given by (3.12), and Qab — C^h '■— |/o^(fl5 ^)| 



Proof. Below we omit 5 for the shortness. Using the relations (5.2 ) and ( |3.7[ ) for -uj 
as in Proposition |5.4| and for the trivial cover, and then (5.5), we get 

Eao...a2„+i [''■(P)] _ F^{aQ, A;a2n+l) ■ FQ{a2n+li<l2n+l) _ 

E+[cr(r)] Fo(ao, A; a2„+i) • F^{a2n+i,a2n+i) 

^2n 



(a2„+i; a2„+i)) ^ 
Z]fe"o(~l)''Pf Go[fc] ■ Fo{ak;a2n+i) ■ (i-Fo(a2n+i; a2„+i))"i 

Pf [ (Gro)j\fe ]o<j"<fe<2n+l 
Efc"o(-l)'PfGo[fc]-(Go)Mn+l " Pf[(Go),,fc]o<j<fc<2„+l ' 
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where and Go are given by (5.4) for corresponding double covers. This can be 
further rewritten as 



ao...a2r.+ lW{^)] 

E+[a(r)] 



Pf 



(Gct 


)j,k (Go) 






[(Go 


)j,k (Go)o,2n+l. 


0<j<k<2n+l 


Pf 


(Go)i,fe 






. (Go)o,2n+l. 


0<j<k<2n+l 



Note that definition (5.4) and convergence (|3.11|) yield 
(Gro)j,fc _ F^{aj;ak)Fo{aj;aj) 



■,7m) 



as 5 



0. Similarly, ([5^ and ( |3.13[ ) give 
(Go)j,fe _ \Fo{aj;ak)\ 



(Go)o,2n+l 

and the factors Caoa2„+i 



|^o(ao; a2r,+i)| |/o^(ao, a2n+i)| 
cancel out in the ratio of Pfaffians. 



□ 



Remark 5.11. (i) One may assume that fQ{a,b) is normalized by (3.8); in this 

case one has a conformal covariance rule Cab — C^{a) <t>(b)\'l^' (^)\^^'^\'^' (^)\^^'^ ■ "^^^ 
linearity of PfafRans implies that the normalization is not important, and that the 
ratio (5.8) is conformally invariant. 

(ii) Clearly, does not change if one adds boundary singletons inside the domain. 
In particular, if = C+ \ {wi, . . . , Ws} and a, 6 G M, then one has (^t — \b — a| 
(with the normalization given by (3.8)). 

(iii) The formula 



|BG93[ equation (17)]. 



was predicted by means of Conformal Field Theory, see 



6. Explicit computations in the half-plane 

In this section we explicitly compute the holomorphic spinors /, 
By (3.12), this immediately gives us the quantities 

and, by conformal invariance of i?'s, all the limits 



C+\{wi,...,Wr,i} 



(6.1) 



lim 



1?((/>(wi), . . .,(j}{Wm)) : 



s'To E+[a(u;f)...a«)] 

for discrete domains {Q^;a^ ,b^) approximating an arbitrary simply connected do- 
main (fi; a, b) and inner faces Wj tending to Wj £ ft, where (f> : C+ is a 
conformal map such that ^(a) = oo and 0(6) = 0. Due to the Pfafhan formula 
(5.8), this result easily extends to the case of 2n+2 marked boundary points. 

It is convenient to choose a = oo. Specializing the boundary value problem 
(a°)-(d°) and (3.8) to the case = C+, a = and using the conformal map 
z —z~^, we get the following conditions for / := /^"''^^"'^' "''"'"^ (oo, •): 



(a°) / is a spinor in C_|- \ {wi, . . . , 

(b°) /(C) e K for any ( e M; 

(c°) ifiz)r^O{\z-w,\-')asz 

(d°) f{z) = 1 + O(z-i) as z^oo. 



j} branching around each of wj; 

Vj and les z=wj{f{z))'^ G iM+ for all j; 
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Note that for m ~ (i.e., the trivial cover), we have an obvious solution /o = 1. 
In order to find /, we introduce an auxiliary spinor 

fioi,....w^{z) ■■= 6^1(2) Sn,^{z) , B^,{z) := ^ ■ 

[[Z — Wj(Z — W)\ 2 

Note that it satisfies (a°), (b°) and (d°). Moreover, it has real zeros at tj = Rewj; 
thus, (a°), (b°) and (d°) also hold for the product fwi,...,w^{z) • g{z), where g{z) is 
any function of the form 

g(z) = l + ^^, AfcGM (6.2) 

(and if some tj coinside, we can add higher-order poles, so that g{z) is always 
a linear combination of m linearly independent functions). We are looking for 
parameters Xj such that fwi,...,w^{z) ■ g{z) satisfy (c°) as well. Denote 

Rk ■■= [-2ireS;,=„,(/^,i,._^„^(z))^]5. 

Then, the condition (c°) for /tui,....™^ (z) • g{z) can be restated as 

™ R 

Im [Rk ■ giwk)] — IniRk + Xj ■ Im =0 for all fc = 1, . . . , m . (6.3) 



Note that (6.3) is an m x m linear system in Xj. We argue that this system is 
always non-degenerate. Indeed, if A° is a solution to the corresponding homoge- 
neous system, then f(z) := EjLi ^^/i^j ~ ^)] ' f'wi,....w^iz) is a spinor satisfying 
(a°)-(c°) and such that f{z) = 0{z^^) as z — > 00. But it follows from the proof of 
Lemma 3.10 that any such spinor is identically zero (after mapping to a bounded 
domain, the condition at infinity yields boundedness near a). Thus, A^ = 0. 

Taking into account that one has explicitly 

Rk = (Imu;,)^ -liB^^iwk) = (Imz«,)^ • 11 . 



and solving the linear system ( 6.3 ) , one obtains Xj as ratios of certain explicit mxm 
determinants (depending of Wj 's) , and then the ratio of spin correlations (6.1) is 
given by 



j = l 3 J j = l 



Re Wj 



Here we have used conventions of Remark 13.71 to determine the sheet of the double 
cover of C+ \ {wi, . . . ,Wn\ (that is, signs of the square roots): in order to get 
the value /(O), one starts with the value f{oo) = +1 and continuously moves the 
boundary point to along the counterclockwise boundary arc (—00, 0), thus arriving 
to Byj.{0) = {—Rewj)/{—\wj\) . The other way to fix the sign is given by taking 
all Wj close to the boundary arc (0; -|-cx)): this should yield a positive correlation 
as that arc carries "4-" boundary conditions. 

Remark 6.1. In particular, for a single point w G C (i.e., for the case m = I) one 
has t ^Rew, R= (Imw)^/^ e M, thus A = and 

Re w 

i!)(w) — - — — — cos [Trhmc, (w,M-)] . 
\w\ 
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Since harmonic measure is conformally invariant, this imphes 
if {V,^; , b^) approximate a, b) and faces tend to an inner point w £ fl. 



cos [7rhmf2(i/7, (ab))] , 
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